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ABSTRACT 

We study the duration and variability of late time X-ray flares following 
gamma-ray bursts (GRBs) observed by the narrow fleld X-ray telescope (XRT) 
aboard the Swift spacecraft. These flares are thought to be indicative of late 
time activity by the central engine that powers the GRB and produced by means 
similar to those which produce the prompt emission. We use a non-parametric 
procedure to study the overall temporal properties of the flares and a structure 
function analysis to look for an evolution of the fundamental variability time- 
scale between the prompt and late time emission. We flnd a strong correlation 
in 28 individual x-ray flares in 18 separate GRBs between the flare duration and 
their time of peak flux since the GRB trigger. We also flnd a qualitative trend 
of decreasing variability as a function of time since trigger, with a characteristic 
minimum variability timescale At/t = 0.1 for most flares. The correlation be- 
tween pulse width and time is consistent with the effects of internal shocks at 
ever increasing collision radii but could also arise from delayed activity by the 
central source. Contemporaneous detections of high energy emission by GLAST 
could test between these two scenarios, as any late time X-ray emission would 
undergo inverse Compton scattering as it passes through the external shock. The 
proflle of this high energy component should depend on the distance between the 
emitting region and the external shock. 

Subject headings: gamma-rays bursts — X-rays: general — high energy: analysis 



1. Introduction 



One of the most unanticipated results to come from the Swift spacecraft (jCehrels et al. 



20041 ) is the wide variety of X-ray behaviors observed in the early afterglows of gamma-ray 
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bursts (GRBs). As of January of 2007, Swift had detected 206 GRBs and had observed a 
subset of > 90% of t hose events with the spacecraft's narrow field X-ray telescope or XRT 
(IBurrows et al.ll2005l ) . Of these events, > 90% show temporal properties that deviate from 
the simple post coohng break powerla w decline that had been seen at late tirn es 3 x 10^ 
seconds) by previous spacecraft (e.g., Frontera et al. 2000; Gendre et al. 20061 ). Afterglows 
with simple powerlaw declines that extend fr om a few ~ 10^ secon ds to several days after 
a burst are seen, for example GRB 061007 (IMundell et al.l l2006l ). but they constitute a 
far minority of the afterglows observed by the XRT. Instead, most after glows show sharp 



drops in the observed flux immediately following the gamma-ray emission (IBarthelmy et al. 



2005al ). lasting anywhere from ~ 10^ to ~ 10^ seconds post trigg er. This is followed by 



a flat tening of the light curve that can last hundreds of seconds (iGranot. Knigl. &: Piran 
20061 ) before eventually transitioning to the late time powerlaw decay previous observed by 
other spacecraft. Most surprisingly, interspersed among these various components of the 
prompt afterglow emission have been the detections of major re-brightening episodes with 
emission flaring in so me cases several hundred times above the declining afterglow emission 
(IBurrows et al.l 120051). In rare cases, t hese flares have actually surpassed the luminosity of 
the original GRB (IBurrows et al.l 120071 ). 



Numerous papers have bee n published discuss i ng a variety of na e chanisms that could 



produce the 1 


ate time flaring ( 
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2006 
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Falcone et 
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2006; 
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2006 
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Proea & Zhane 
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Lazzati & Pernall2007 


Lee et al. 


2007; 


Lvutiko\ 


^ I2OO6: 


Fan & We 


i 2005). Most of these mechanisms place tight constraints 



on th e timescales on which the their emission can be produced (lloka. Kobayashi fc Zhang 
2OO5I ). The simplest explanation would be that the forward shock powering the afterglow runs 
into ambient density fluctuations as it moves into the surrounding medium (jWang fc Loeb 
2OOOI ). This external shock interpretation has diffi culties explaining th e degree of variability 
that is clearly seen in many of these flares (e.g., IBurrows et al.l l2005l . and below). Simple 
kinematic arguments show that fluctuations due to turbulence of the interstellar medium 
or variable winds from the pr ogenitor are expected to produc e broad and smooth rise and 
decay proflles, with At/t ~ 1 (lloka. Kobayashi &: Zhangll2005l ). Here t is the time since the 
gamma- ray trigger and At is the variability timescale. Sh ocks internal to the relativistic 
outflow (iRees &: Meszaroslll994l : iNarayan fc Paczynskilll992l ). similar to the shocks believed 
to produced the prompt gamma-ray emission, do not suffer from these same constraints 
and could in theory produce variability on much shorter timescales. In the internal shock 
scenario the rise time of an individual pulse is governed by the time it takes for the reverse 
shock to propagate back through the shell. The decay time is largely set by the relativistic 
kinematics, or curvature effects, in which the arrival of off axis emission from a relativistically 
expanding shell is delayed and affected by a varying Doppler boost. 
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Another clue that the flares are produced in a region distinct from the external shock is 
that the temporal decay of the afterglow emission appears largely unaffected by the presence 
of flaring. The temporal index of the afterglow after the flaring activity is typically consistent 
with the pre-flare decay index. Although most bright flaring occurs within one hour of the 
GRB, flares have been observed during of each of the light curve phases described above. If, 
for example, the flare represented the onset of forward or reverse shock emission of a slow 
shell catching up and colliding with the external shock, then these flar es would be expected 
to occur only before the flat energy injection phase. Further more. (IBurrows et al. l2007l) 
points out flaring in one example of a possible "naked burst" ( Kumar &: Panaitescu 200ol ). 
an event which decays rapidly in time and therefore exhibits no evidence of external shock 
emission. This supports the argument that whatever is powering the afterglow is most likely 
not creating the X-ray flares, leaving internal shocks or direct central engine activity as likely 
methods for their production. 

Further evidence that late time X-ray flares might be associated with internal shocks 
comes from their spectral characteristics. First, most of the flares ar e muc h harder than 
the underlying afterglow emission and, as reported by iBurrows et al.l (120071 ). the spectral 
characteristics of the afterglow emission appears unaffected by the fla ring activity, possi- 
bly iii dicating two distinct emitting regions. Second, spectral fltting by iButler &: Kocevski 
(120071 ) has shown that many flares can be well flt by the Band model (iBand et al.l Il993l ) 
that so effectively describes the prompt emission which is largely believed to be the result of 
int ernal shock collision . Furthermore, detailed time resolved spectral fltting of bright flares 



by iButler &: Kocevskil (120071 ) has shown that the spectral break energy Ep^ of the Band 



model, which represents the energy at which most of the photons are emitted, evolves to 



lower ene rgy during t 



le flar e in a way that is very similar to what is seen in the prompt 



emission (Norris et al 



19861 ). The evolution also follows the hardness- intensity correlation 



(IGolenetski et al.lll983l ). a well known relationship observed in the prompt emission that can 
be attributed to the relativisti c effects that produce the decay proflle of individual pulses 
jKocevski. Rvde fc Liang! I2OO3I ). 



If the energy released by this activity is converted to radiation through late time internal 
shocks, then the question remains as to the characteristic radius that these internal shocks are 
occurring as well as the delay in their ejection. Either the central engine is still functioning 
and emitting shells at very late times, or the flnal few shells of the original outflow, which 
were emitted along with the shells that created the prompt emission, catch up with each other 
only after a long delay due to a small relative difference in their bulk Lorentz factor F. The 
flrst scenario could essentially produce shell collisions at any radius, as the delayed arrival 
of the flares would, in th i s case , primarily reflect the time that the engine was dormant 
(iKobayashi. Piran. fc Saril 119971 ). The second scenario, predicts that the late time flares 
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should occur at a radius that is significantly larger than the radius at which the prompt 
emission was created, with their delayed arrival being a result of the shells' time of flight 
before colliding. This second scenario leads to a very specific and testable prediction, namely 
that the width of individual pulses of emission should beco me broader and less variabl e when 
originating from shells of increasing collision radii Re- iRamirez-Ruiz fc Fenimord (120001 ) 
tested for this pulse width evolution in the light curve profiles of BATSE events and found 
no evidence for any such effect. They concluded that the prompt emission observed by 
BATSE must have been produced over a small range of Rc from the central engine and 
that no significant deceleration of F could have occurred over the duration of the observed 
activity. 

The goal of this paper is to extend the gamma-ray pulse width analysis to the late 
time flaring X-ray emission following GRBs. The public catalog of Swift XRT flares (see 
also Chincarini et al. 2007) represents the first dataset to test the internal vs. external 
shock scenario for this flaring activity. Whereas previous studies were limited to prompt 
emission occurring less than 100 seconds after trigger, the late time X-ray flares give us 
the opportunity to test for pulse width and variability evolution out to, in some case, 1000 
seconds after the trigger of the GRB where this effect may be more pronounced. We provide 
a simple derivation of the expected pulse width evolution in both small AF and delayed 
engine activity scenarios in §2, followed by a discussion of our data reduction techniques in 
§3 and results in §4. We find evidence for pulse width evolution in 28 flares as well as a 
qualitative trend of decreasing variability as a function of the flare's time of peak flux. We 
discuss the implicati ons of our observations in ^5. This work expa nds upon and formalizes 
our previous reports (IKocevski. Butler fc Bloomll2006l : lButleiil2007l ) of the discovery of pulse 
width evolution. 



Pulse Width Evolution 



The standard fireball model postulate s the release of a larg e amount of energy by a 
central engine into a concentrated volume (ICavallo fc Reeslll978r). which causes the result- 
ing outflow to expand and quickly be come relativistic (jPaczynskil 1 19861 ). In the internal 
shock scenario (IRees fc Meszarod Il994l ). this outflow is assumed to be variable, consisting 
of multiple shells of differing bulk Lorentz factors F. These shells propagate and expand 
adiabatically until a faster shell collides with a slower one, causing the shells to coalesce and 
convert a significant fraction of their kinetic energy into radiation, most probably through 
optically thin synchrotron radiation. The resulting pulse profile that is observed is a con- 
volution of two distinct timescales. The rise time of the pulse is largely due to the time it 
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takes for the reverse shock that is induced by the colhsion to cross the width of the faster 
shelL The decay time, on the other hand, is governed mainly by angular and kinematic 
effects where off axis emission is delayed and effect ed by a varying Doppler boost due to 
the curvature of the relativistic shell (see Figure 1 in iKocevski. Ryde &: Liang ISoJ). As 



a 

result, the decay time can be, and in most cases is, much longer than the rise time, leading 
to an asymmetric pulse profile. The combination of these two timescales (the shell crossing 
time and the angular time) naturally explains the so called "fast rise exponential decay" or 
FRED pulses that are so ubiquitous in prompt GRB emissionjl] 

If we examine these two timescales in more detail, we can see that the rise time is 
primarily a thickness effect and can be expressed as Atrise = 5R/c{[32 — Prs), where 6R and 
(32 are the thickness and velocity of the second shell that is catching up to the first and jSrs 
is the velocity of the reverse shock. If both the slow and fast shells have Lorentz factors of 
roughly the same order ~ F, then the resulting rise time is of order ~ 6R/c. Because the 
merging shells are traveling forward at a velocity very close to the speed of light (F ^ 1), the 
resulting coalesced shell keeps up with the photons that it emits. Therefore, any emission 
activity over a fixed duration will appear to an outside observer to be compressed in time 
by a factor of 1/2F^, where F^ is the resulting Lorentz factor of the merged shell. The 
observed rise time can therefore be written as 

At 



2cF 



So given a sufficiently large F^, internal shocks can essentially produce variability along the 
line of sight on arbitrarily short timescales. 

Angular (or curvature) effects have the opposite effect, causing a broadening of the 
overall emission profile that can quickly come to dominate the observed pulse shape. The 
decay timescale is essentially the difference in light-travel time between photons emitted 
along the line of sight and photons emitted at an angle 6 along a shell of radius R. This can 
be stated as 

i?(l -cosA^) _ R{A9Y _ R 
~c ~ Yc ~ 

Where the last step assumes that the shell is moving with sufficient velocity such that the 
solid angle accessible to the observer is limited by relativistic beaming and thus given by 
AO ~ 1/F^. Therefore comparing Equation [1] and Equation [21 we can see that curvature 



At. = l^l^i^^^ ^ 1^ ^ ^ (2) 



^Here we assume that the intrinsic coohng time Ate of the sheU is insignificant compared to the duration 
of the shell crossing /Sir and angular /Sid timescales because of the magnetic field strength required to 
produce the gamma-ray emission 
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effects become important whenever the radius of the shell exceeds the shell thickness, which 
is true for all but the earliest moments of the shell's expansion. 

The significance of Equation [2] is that angular effects should scale linearly with the radius 
of the emitting shell and therefore pulse durations should become broader as shell collisions 
occur further from the central engine. If the flares are the result of multiple shells that have 
been ejected almost instantaneously (or at least within a timescale that is small compared 
to the overall GRB duration) but collide at very late times due to a small dispersion in 
Lorentz factors, then one would expect that these late collisions would occur at greater 
radii. In this scenario, we can replace the radius of the shell in Equation [2] with the time 
t since the ejection of the first shell by noting that the observed radius of a spherical shell 
expanding with f ~ c can be approximated as i? ~ ctT"^, where the extra factor of is due 
to relativistic corrections, leading to 

At. ^ I (3) 

Therefore, the late shock scenario would predict a linear correlation between a shell's time of 
flight and the resulting pulse duration, independent of the Lorentz factor of the shell. This 
relationship between th e pulse duration the time sinc e the ejection of the internal shocks 



has been noted before. iFenimore. Madras. Nayakshiru (119961 ) found, through a much more 



detailed derivation, that a pulse's FWHM should scale ro ughly as 0.26Tn to 0.19Tn as the 



low energy power law index a varies from 1 to 2. Similarly, lloka. Kobayashi fc Zhang] (120051 ) 
derive that the variability of flares that result from refreshed shocks should be limited by 
At ^ tp/4. In each case, flares occurring at larger radii are expected to produce broader 
pulse durations. 

This relationship between Atd and t is modified if there is an intrinsic delay Aten^me in 
the ejection of the subsequent shells by the central engine. If we imagine two shells emitting 
at time zero and time Atengine, provided the Lorentz factor of the second shell r2 > Fi, the 
Lorentz factor of the first shell, the shells will collide at time 



F^Af 

, _ ^ 1 ■^''engine / ,n 

If the the shells have equal mass, which corresponds to the maximal efficiency for conversion 
of kinetic energy into radiation, energy and momentum conservation lead to a merged shell 
with Lorentz factor F^ = y/TiT^. The timescale over which the shell emits will be governed 
by the longest timescale of At^, At^, or At^ the angular, the radial, or the cooling timescale, 
respectively. The angular and radial timescales can both be given as: 



The time at which a flare is observed will be tc, and the observed duration will be At = 
tcT i/r2 ~ tc/2 for an efficient collision with r2 = 2ri. For this r2/ri ratio, the flare duration 
At is related to the duration at the central engine by At = Atengme/6 ~ Atcngine- Therefore, 
if there is any appreciable delay in the ejection of relativistic material from the central engine, 
the resulting pulse shape will not necessarily reflect the shell radius, but rather the intrinsic 
delay between the ejection of the two shells. Any correlation between pulse shape and time 
of peak flux must then be attributed to the activity of the central engine. 



3. Data & Analysis 

We select a subsample of 28 bright 10 cts/s) flares that are fully time-sampled 
(i.e., no gaps in their light curves) in 18 separate GRB afterglows observed by XRT. The 
Burst Alert Telescope (BAT) and XRT data were downloaded from the Swift Archiv^ and 
processed with version 0.10.3 of the xrtpipeline reduction script and other tools from the 
HEAsoft G.O.cjfl software release. We employ the latest (2006-12-19) calibration files available 
to us at the time of writing. 

The reduction from cleaned event lists output by the xrtpipeline code and from 
the H EAsoft BAT softwa r e to science ready light curves and spectra is described in de- 



tail in iButler fc Kocevskil (120071 ). The bright XRT flare data are taken overwhelmingly in 
windowed-timing (WT) mode, which mandates special attention to bad detector columns. 
As the spacecraft moves, a significant and time varying fraction of the source flux can be los t 



if source counts fall on the bad columns. To account for this (see, iButler fc Kocevskil 120071 ). 
we calculate exposure maps for the WT mode on a frame-by-frame basis. We accumulate 
0.3-10.0 keV counts in each light curve bin until a fixed signal-to-noise [S/N) of 3 is achieved. 
All of the resulting light curves to which we apply our analysis are publicly availabl^. A 
composite light curve plot showing all 18 GRBs in our data set is shown in Figure 1. 



^ftp: / /legacy.gsfc. nasa.gov/swift /data 

■^http: / /heasarc. gsfc.nasa.gov/docs/software/lheasoft/ 

^http: / / astro.berkeley.edu/~nat/swift 
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3.1. Flcire Duration Measures 

The first step in our analysis consists of measuring the global fiare duration timescales. 
Because we have found no one functional form (e.g., Gaussians) to adequately fit the X-ray 
flare time profiles (Figure 1), we employ non-parametric duration estimators. We consider 
the flare Tgo duration as the time required to accumulate between 5% and 95% of the flare 
counts. We also deflne a rise time as the time between 5% accumulation and the time of count 
rate peak. Errors on these quantities are determined from the non-parametric bootstrap (i.e., 
by recalculating the quantities for data simulated using the measured data and errors). 

A bias affecting these duration measures (and probably all durations measures) is the 
unknown background under the flares. As discussed above, studies have shown that a power- 
law decaying background likely does exist. However, it cannot cleanly be measured in many 
of our events and not at all for events which suffer from data gaps. In an effort to avoid 
such biases in our duration measurements, we have restricted our analysis to flares that are 
typically 2—3 orders of magnitude above background. 

3.2. Flare Variability Measures 

The flare duration and the component rise and decay times are gross measures of vari- 
ability. In addition to this information, we attempt to measure the flner timescale fluctua- 
tions in the light curves which may prove important for inferring the size and nature of the 
flare's emitting region. Several methods for measuring signal power versus time scale have 
been applied to astronomical inquires and in GRB research in particular. Several authors 
(e.g.. Belli 1992; Giblin, Kouveliotou, & van Paradijs 1998; Beloborodov, Stern, & Svensson 
2000) employ the Fourier power spectral density (PSD) to study time variations in GRB 
light curves. The autocorrelation function (ACF), which is simply the Fourier transform of 
this PSD, has been used to demonstrate a narrowing of GRB pulses with increasing energy 
band (e.g., Fenimore et al. 1995). Below, we utilize the first order structure function, which 
is directly proportional to the ACF and has a rich heritage in the study of quasar time 
histories (e. g., Simonetti, et al. 1985; Hughes, AUer, & AUer 1992). 

Because the light curves of flaring sources are by deflnition non-stationary signals (i.e. 

signals whose frequency content changes with time) which exhibit sharp discontinuities, 
Fourier transforms do a particularly poor job of accurately measuring their power on both 
short and long timescales. Furthermore, they offer no ability to distinguish the temporal 
variations of specific spectral components (i.e. the time at which a characteristic frequency 
changes in a light curve) They are also somewhat more prone than ACF methods to aliasing 
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effects due to irregularly time-sampled data. 

Instead of constructing a PSD using superpositions of sines and cosines, we can perform 
the equivalent analysis by constructing a scaleogram through the use of a discrete Haar 
wavelet transform. The Haar wavelet is the simplest possible wavelet, consisting of a step 
function, and has been previously exploited to "denoise" GRB light curves (e.g., Kolaczyk 
& Dix on 2000) and to infer milli-second va riability during the first seconds of bright BATSE 



GRB flWalker. Schaefer k Fenimorell2000l ). 



As described in more detail in the Appendix, we calculate the structure function from 
Haar wavelet coefficients as: 



2 



t/2Ai-l 

At/t J2 X2iAt?. (6) 

i=0 

where Xi is the natural logarithm of the observed XRT count rate in bin i at time t, and 
At is the timescale (or time "lag" ) between successive bins. The bar over the denotes an 
averaging with respect to shorter timescales, which is accomplished by the discrete wavelet 
transform (see, e.g.. Press et al. 1992). If this averaging were not performed, Cx At would 
be equal to the structure function SF =< (Xj+At — Xi)'^ >= 1 — ACF. Instead, we have an 
estimator for SF, which ends up being far easier to interpret, as we discuss in the Appendix. 



4. Results 

4.1. Pulse Broadening in an Individual Event 

A composite BAT and XRT light curve for GRB 060714 is shown in Figure 2. The 
red solid line represents a multiply-broken powerlaw fit to the light curve. The inflection 
points in the fit allow us to measure the boundaries and durations of the individual pulses 
within the signal. Each pulse is delineated by the short-dotted lines with the corresponding 
pulse duration labeled below the light curve. A general trend can be seen in which the pulse 
durations become broader as the burst progresses, with the shortest activity occurring early 
in the event. 

The bottom panel plots the minimum timescale At for which ax,At is at least 3a above 
the fioor expected from Poissonian fiuctuations. Consistent with the trend seen in pulse 
duration, this minimum variability timescale, which is calculated without fitting the data, 
increases roughly as a powerlaw as the burst moves from early gamma-ray emission to late 
X-ray emission. To show that this increase in the variability timescale is not simply due to 
an increase in the data binning as the burst fades, we have plotted the time binning as a 
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short dotted line in the bottom panel. For most of the event [t < 200 seconds), the timescale 
below which little or no significant power exists is at least an order of magnitude higher than 
the resolution of the light curve allowed by the binning of the signal. Overall, the bottom 
panel shows that the very fast time variability associated with the prompt GRB emission 
dies out at late times. 

Although the general broadening of the pulse durations seen in Figure 2 can be detected 
within the separate BAT and XRT light curves, the comparison of the pulse durations can 
only be qualitative when considering a light c urve that spans both detectors. This is because 
GRBs are typically wider at lower energies (IFenimore. Madras. NayakshinI Il996l ) , a direct 
result of the evolu t ion of their spectral break energy Epk to lower energies. We shown in 
Butler &: Kocevskil (120071 ) that the X-ray flares typically have Ep^ in the X-ray band, while 
the earlier GRB emission has Ep^ in the gamma-ray band. Therefore pulses are expected to 
be intrinsically broader in the 0.3—10.0 keV bandpass of the XRT than the higher 10—100 
keV observed by BAT. 



4.2. Pulse Broadening in the Sample Taken as a Whole 

To eliminate the pulse broadening between separate energy bands, we limit our quanti- 
tative comparison of pulse durations (both within a single GRB and across our entire sample) 
to measurements made using only the XRT data on each event. This comparison is shown in 
Figure 3, where we plot pulse duration versus time of peak flux for our entire sample of XRT 
observed GRBs. The flares associated with each GRB are represented by the same color and 
symbol, with several GRBs exhibiting multiple flares throughout their early afterglow. 

As a whole, the sample shows a clear correlation between the flare duration and the time 
of peak flux since the GRB trigger. The resulting correlation strength is Kendall's tk = 0.7, 
with a signiflcance of 10"''. The slope is consistent with linear, implying At oc tp, which 
cancels out the effects of cosmo logical redshift. We flnd no signiflcant correlation between 
duration and redshift {tk = 0.2, signif.= 0.2), further ruling out cosmological time dilation 
as the source of this correlation. 

Roughly, half of the events with multiple flares (those plotted in color in Figure 3) 
show a trend toward increasing duration with observation time. The other half show an 
anticorrelation. The pulse durations, time of peak flux, and rise times for al the flares in our 
sample can be found in Table 1. 
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4.3. Haar Structure Function View 

Figure 4 shows oxAi versus At and ax At versus At/t for the ensemble of flares under 
stud}0. In this scaleogram plot, we show only Sex excesses over the power associated with 
Poisson fluctuations and report lower values as Scr upper limits. 

An X-ray flare is an emission episode uncorrelated in time to the afterglow flux prior to 
and after the flare. During the flare and on timescales short relative to the flare duration, 
the flux will be highly correlated in time and there will be a linear rise in ax, At- This can be 
observed to arbitrarily short timescale if the fading powerlaw tail of a flare is measured with 
very high S/N. On the other hand, as we describe in more detail in the appendix, correlated 
behavior in the light curve flattens the structure function, and this provides a direct measure 
of the flare timescales. 

Consistent with the pulse duration correlation seen in Figure 3, the scaleogram plots 
show a range of important flaring timescale (it = 30 — 300s, which becomes much tighter in 
units dt/t = 0.1 — 0.5. We observe a minimum characteristic timescale dt/t = 0.1. 



The fractional flux variation levels at the minimum timescale are large {ax,At ^ 80%), 
suggesting that the variations correspond to gross features in the light curve. Consistent with 
this interpretation, we observe the flare rise times to have Atrise/t = 0.1 on average (Figure 
8]), and it is likely the sharp flare rises which produce the shortest timescales reflected in the 
structure function turnover. From the linear ax,Ati we can rule out significant flickering on 
timescales shorter than dt/t = 0.1 (or dt = 30s) at very small ^ 3% fractional flux levels 
(Figure 3). We discuss the flare n oise properties as a function of timescale in more detail in 
Butler. Kocevski Sz Bloom J2007h . 



For observation times in the 100 to 1000 second range, dt/t = 0.1 implies emission radii 
Rc ~ 10^^ cm— 10^^ cm, for a bulk Lorentz factor F = 100 (Equation [1]). The observable 
emission is restricted to an angle ~ l/F, implying an effective emitting region of size SR ~ 
i?c/r ~ 10^^ cm — 10^^ cm, c ompared to the typical external shock values of 10^^ cm in the 



first hour or so (jPirarull999l ) 



^We reserve a more det ailed study of the Haar structur e functions of individual flare events and Swift 
GRBs in a separate paper (jButler. Kocevski fc Bloomll2007( ). 
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5. Discussion 

The results from the temporal analysis outlined above provide substantial evidence that 
both the pulse duration and pulse variability of late time X-ray flares evolve with time. 
Both the pulse duration and variability timescales appear to have a narrow intrinsic range 



i n At I t = 0.3 ± 0.2, cons istent with a nar r ow ra nge found independently by iBurrows et al. 



(120071 ) for Atrisc and by IChincarini et all (120071 ) for Tgo. GRB 060714 provides the best 



example of this behavior in an individual event. Several other individual GRBs display a 
similar increasing pulse duration trend among their associated flares, although several bursts 
do not (e.g., GRB 060210). For the bursts with multiple flares, only half show increasing 
flare durations. Each burst event typically shows only 1, sometimes 2 (and 3 in one case) 
separate flares. These multiple flares within individual GRBs are only weakly separated in 
logarithmic time, and hence probe a small range of Rc or tengine which may not allow for a 
clean measurement of time evolution in individual events. 

The linear relationship between At and tp is consistent with the pulse width evolution 
that is expected from the angular effects of late internal shocks at large radii as outlined in ^ 
Because we do not see a significant alteration of the afterglow light curve after the occurrence 
of an X-ray flare, the standard refreshed shock model, in which the trailing shells catch up 
to the leading shell only after the leading shell decelerates due to an external medium, is 
disfavored. Although such a scenario is expected to pr oduce a correlation between the pulse 



width and time of peak flux on the order of At ^ tp/4 (lloka. Kobayashi fc Zhangll2005l ). the 
trailing shells should have the effect of increasing the overall afterglow energy and thus have 
a discernible effect on the afterglow light curve, which is not seen. Therefore, the internal 
shocks producing the flares would have to be occurring behind the leading shock that has 
begun powering the afterglow, with their late occurrence, in this scenario, being due to a 
small relative Lorentz factor between the two inner shells. 

The primary difficulty with this interpretation is the high flux ratio between the prompt 
and late emission, given the relat i vely s mall AF needed to explain the late collision time. As 



shown in detail by iKrimm et al.l (120071 ). the efficiency e of an internal shock in converting a 
system's kinetic energy into radiation scales roughly as e ~ AF^, so the observed flux drops 
quickly as the contrast between the Lorentz factor of the shocks decreases. This posses a 
problem for the flares observed by Swift, as many exhibit peak fluxes that are significant 
fractions of, and in some cases comparable to, their associated prompt emission. The small 
AF scenario would require an extremely large total amount of kinetic energy to remain in the 
system after the release of the prompt emission, given the low efficiency of the late collisions. 
These late and highly energetic shocks would, after producing the flaring activity, eventually 
collide with the external shocks and affect the observed afterglow light curve, something that 



- 13 - 



is not seen in all events with flares. 



Alternatively, if the late nature of the X-ray flares is due to a significant delay in the 
ejection of late shells by the central engine, then the necessity of a small AF is eliminated, 
alleviating this efficiency constraint. As described in ^ the arrival time tc and pulse width 
At would then directly refiect the activity of the central engine. Therefore, in this scenario, 
the correlation between tc and At would require an explanation intrinsic to the powering 
and/or reactivation of the central engine at late times. Several authors have suggested 
mechanisms by which the central engine could be active at late times, most involving late- 
time fallbac k material or a lon g lived accretion disks around a central black hole. A model 
proposed by lKing et al.l (120051 ) suggests that the late-time activity could be attributed to the 
fragmentation and accretion of a collapsed stellar core resulting in a sporadic release of e nergy 
rather than the classic view of a single cataclysmic event. Similarly. iPerna et al.l (120061 ) have 
proposed a viscous disk model in which the late-time activity is due to re-energization by 
material that falls in from a range of initial radii toward the accreting black hole. In this 
scenario, the correlation between tc and At would be due to the range of radii from which the 
accreting material was falling. Material at large radii, if continuously distributed throughout 
its orbit, would take longer to fall back onto the central black hole and would do so over a 
longer duration, due to its larger orbital circumference. 

Th ese models are no t without their own share of difficulties. The simple fragmentation 
models (IKing et al.ll2005[) are inco nsistent with the implication of the spectral evolution seen 
in may flares (IKrimm et al.l 120071 ). Similarly, the viscous disk model requires a continuous 
distribution of material at discrete orbits to account for the episodic nature of the flares as 
well as an extremely long lived, and hence low viscosity, accretion disk to explain flares at 
1000 seconds after the original collapse. 

It cannot be completely ruled out that the observed time ev olution is due to spectral 
evolution or the sup erposition of multiple flares. Consider 060124 (IButler &: Kocevskil 120071 : 
Romano et al.ll2006l ). in which the flares may in fact be the prompt emission, because the 
faint BAT trigger may be a pre-cursor. At high energies, the flrst XRT "flare" resolves into 
2—3 shorter timescale BAT flares, which are blurred together in the XRT. We note that a 
shift of time origin for 060124 from t ~ s to t ~ 300 s, corresponding to a shift in origin 
from the pre-cursor to the flare start, does not lead to a violation of the At and tp correlation. 
Although, if we used the BAT flare durations, the correlation could be violated. 

This indicates that spectral considerations are important, and that we are likely measur- 
ing in the XRT (in some cases) a pulse superposition. The duration which increases in time 
still appears to measure the duration of major emission activity, however, it is not clear that 
these are individual pulses. We know that spectra of late time flares are evolving strongly 
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( iButler &: Kocevskil 120071 ) during the flares. However, we observe only a weak correlation 
between peak time and hardness, indicating that there is a diversity of flare spectra at each 
epoch. 

Another important concern involves the powerlaw background onto which most of these 
flares are superimposed. Although we have not attempted to subtract the background from 
the events in our sample (because the backgrounds are not well defined), this should not 
dominate the observed correlation. We have selected the brightest flares for analysis, which 
have peak fluxes orders of magnitude greater than the underly ing background flux. The 



correlation is also strong for measures of duration like T50 or the lReichart et al.l (l200ll ) T45, 
which are largely insensitive to pulse tails. Finally, we note that the flare rise time also 
strongly correlates with the peak time tp as shown in Figure 5. 

Barring any of these selection and/or analysis effects and assuming that the pulse width 
evolution is real, one possible test to distinguish between the late internal shocks with small 
contrasts AF and direct central engine activity may come from contemporaneous high en- 
ergy emission during the X-ray flares. If the internal shocks creating the flares are occurring 
behind the external shock, then one would expect the X-ray photons to be boosted to higher 
energies by a factor of F fq through inverse C ompton (IC) scattering as they pass through 



the external shock flRybicki &: Lightmanlll979l ). The soft X-ray 10 KeV photons associated 



with the X-ray flares could easily be boosted into the 1—100 MeV range depending on the 
Lorentz factor of the external shock. The temporal profile of this high energy component 
should dep end heavily on the distance be hind the external shock at which this emission 



originated (jWang. Zhuo. fc Meszarosll2000l ). as the duration of the IC component will reflect 
the geometry of the external shock, roughly i?/2F^c. The ratio between the flare duration 
and the IC component's duration should approach 1:1 as the radius of the internal shock 
producing the flare approaches the external shock radius. Internal shocks that result from 
delayed central engine activity do not necessarily have to be at large radii to produce the 
longer observed durations. Therefore, larger IC component to flare duration ratios are ex- 
pected for flares produced from small radii collisions. Even if these late-time collisions at 
small radii have intrinsically longer durations, as suggested by late central engine activity 
models, the additional light travel time from the origin of the late time flares to the external 
shock as it expands may make this change in duration ratios measurable. Such a test for 
contemporaneous high energy emission will be aptly suited for the upcoming GLAST mission 
which will be sensitive to photons up to > 300 GeV. 
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7. Appendix 

We describe here the mathematical representation of a Haar wavelet and its use in the 
construction of a scaleogram closely related the the ACF and first order structure function 
SF. 

Given T successive data bins Xi, we define the Haar wavelet coefficients hi^i on scale 
At = 1 as 

^ X2i+, - X2i, i = 0,...,T/2-l. (7) 
At the same time, we can calculate the signal smoothed over a 2 bin scale At = 2: 

Xi,2 = liX2i+i + X2i), i = 0, T/2 - 1. (8) 

By successively differencing and smoothing the signal on dyadic scales At — 1,2, A, etc., we 
build up the discrete Haar transform (see, also. Press et al. 1992): 

Km = X2i+iAt - X2i,At: i = 0, ...,T/2At - 1. (9) 

If the Xi are uncorrelated with equal variance, then the /ij^At will be approximately linearly 
independent. We form a Haar scaleogram by averaging the Km a-t each scale At: 

t/2M-l t/2M-l 

(T%M-^t/t J2 K^t^At/t iX2i+l,M-X2i,Mf. (10) 

1=0 i=0 

In practice, we calculate this average as an average weighted by the data measurement 
uncertainties, Wi — l/crj),i- 

This quantity, also known as the Allan (1966) variance, is closely related to the structure 
function SF —< {Xi+M — >, where < ... > denotes an average over the data. Unlike 
<^l:,At the quantity SF is calculated without averaging the data on scale At before differencing 
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on that scale. This leads to a scaleogram with correlations (even for uncorrelated input data) 
between nearby data bins. The uncorrelated scaleogram aj^ therefore easier to fit and 
interpret, while both scaleograms have similar shapes for a wide variety of noise models. 

7.1. Flare Ensemble Haar Structure Function 

Because the Haar wavelets encode signal scale information as a function of time, it is 
possible to calculate c^x At ^"^^ arbitrary time sections of a light curve (e.g.. Figure 2) or for 
the full light curve. 

To make useful scaleogram plots for multiple GRB flares (e.g.. Figures 4), we place the 
times series data end-to-end and perform the Haar transform as though the data were binned 
on an even time grid. Transform coefficients formed by differencing data from separate events 
are discarded. By saving the actual time since GRB trigger t and time bin width At for each 
wavelet coefficient, we can then rebin the coefficients in time on a dyadic grid starting with 
the minimum bin size. In this fashion, it is possible to plot statistically independent ax,At 
points versus the physically meaningful At or At/t. 

For Xi in Equations (THHl we use the natural logarithm of the XRT count rate. Because 
the counts have been binned to a fixed S/N ratio, the error in Xi is approximately constant 
((Td ~ 1/3). The natural logarithm is also useful because powerlaw flux variations lead to 
a "zero-flaring" scaleogram with ax,At oc At, as can be seen from a Taylor expansion of 
the flux in time. Also because we are working with the logarithm of the count rate, ax, At 
can be interpreted as a root-mean-square (RMS) fractional variation in the flux F (i.e., 
6X ^ 6F/F). 

7.2. Structure Function Interpretation 

Following the discussion in Hughes, AUer, & AUer (1992), on short timescales, the 
scaleogram ax,At asymptotes to a^, where is the data measurement uncertainty. Because 
we know a^, we can subtract this flattening out. (This is typically not possible for SF due 
to the introduction of correlations in the data.) From the Cauchy-Schwarz inequality!^, the 
scaleogram increases with increasing time lag. It eventually saturates to a characteristic 
signal level o"signai at time ^ Tgo, once we begin to run out of correlated variations in the 



^Recall that the Cauchy-Schwarz inequahty states that \{x,y)\'^ ^ {x, x) ■ {y, y) and that the two sides are 
equal only if x and y are Hnearly dependent. 
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signal. 

On intermediate timescales, the slope of crx,At depends on the shape of the hght curve 
and on the noise spectrum of possible low-level or unresolved flares. If the light curve is 
correlated on these timescales, which is to say smooth on these timescales, ax,At will increase 
as At. If, however, the light curve is dominated by the sum of slowly decaying responses to 
low level flares, a characteristic "flicker noise" spectrum {PSD{f) oc 1//) may result and 
o'x,At ^ ^^"^ Hence, we can test for flaring as a function of timescale by measuring powerlaw 
o'xAt slopes less than unity. 

The fading powerlaw tail of a flare measured with infinite S/N would produce a statis- 
tically significant (Tx,At for arbitrarily small At. These timescales, where crx,At oc At, are 
therefore uninteresting. However, the beginning of a ax,At At° phase yields a physically 
meaningful timescale for the flaring. The breadth of this phase indicates the range of At 
present in the light curve. 
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8. Figure Captions 

Fig. 1. - The XRT count rate (cts/s) plotted vs. time since trigger for all 28 flares in 
18 separate GRBs. The light curves in this plot are rebinned to S/N — 10. A qualitative 
trend between pulse width and time of peak flux can be seen by inspection. 

Fig. 2. - Top Panel. A composite BAT and XRT light curve for GRB 060714 showing 
an increasing pulse duration as a function of time. Bottom Panel. The minimum variability 
timescale in the composite light curve (with power that is at least 3a above that which is 
expected from Poissonian fluctuations). The variability of the light curve increases with 
time, roughly as a powerlaw of ATmin oc T^-^^^-^, 

Fig. 3. - The pulse duration Tgo versus time of peak flux Tp for our entire sample 
of XRT observed flares. Multiple flares from individual GRBs are displayed with a unique 
color-symbol combination, whereas GRBs with only one flare are represented by a black 
diamond. A strong trend {taux = 0.7) between pulse width and the time since trigger, as 
measured in the observer frame, is clear from the data. Only half of the GRBs with multiple 
flares display a similar increasing pulse duration trend between their associated flares. We 
conclude that the observed pulse width evolution only becomes apparent when examining 
durations that cover a broad temporal range. 

Fig. 4. - Haar wavelet scaleogram ax,M versus timescale At (Panel A) and At/t 
(Panel B) for the ensemble of flares under study. The expected level for Poisson noise 
has been subtracted out. Because (Tx,At is calculated from the natural logarithm of the 
XRT count rate, it can be interpreted as a measure of RMS fractional flux variation versus 
timescale. The scaleograms reach maximum and turn over on timescales At f=:i 30 — 300s 
and At/t k, 0.1 — 0.5, indicating that the flaring occurs on these characteristic timescales. 
Signiflcant (> 3-sigma level) variability is observed on timescales At ^ 3s and At/t ^ 0.01, 
however, ax, At oc At (dotted red curves) indicates that this variation is due to flaring on 
intrinsically longer timescales. 

Fig. 5. - The flare rise time T,. plotted vs. the time of peak flux Tp. As in Figure 3, 
multiple flares from individual GRBs are displayed with a unique color-symbol combination, 
whereas GRBs with only one flare are represented by a black diamond. An increasing trend 
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similar to that seen between T90 and Tp is evident in the data. The observed rise times are 
largely insensitive to the effects of background subtraction. 
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Table 1. Pulse Duration & Time of Peak Flux For Our Entire Sample 



GRB 


Time Region 


^90 






T- 

rise 




(s) 




(s) 




(s) 


(s) 




050502B 


400.0 - 


1200.0 


358.1 ± 3.2 


784.6 ± 23.7 


203.5 ± 23.7 


050607 


250.0 - 


- 600.0 


165.5 ± 18.0 


312.6 ± 3.2 


28.8 ± 9.6 


050713A 


95.0 - 


- 150.0 


39.3 ± 


0.5 


122.8 ± 5.3 


20.7 ± 5.3 


060111A 


200.0 - 


- 500.0 


191.3 ± 


2.5 


279.6 ± 15.4 


43.4 ± 15.4 


060312 


100.0 - 


- 200.0 


56.0 ± 


2.5 


113.2 ± 1.7 


11.4 ± 1.7 


060526 


230.0 - 


- 450.0 


128.9 ± 


0.9 


251.4 ± 2.8 


4.8 ± 2.8 


060604 


120.0 - 


- 200.0 


63.9 ± 


0.5 


136.8 ± 3.4 


11.0 ± 3.5 


060904B 


140.0 - 


- 300.0 


105.4 ± 


0.8 


182.3 ± 6.6 


30.5 ± 6.7 


060929 


470.0 - 


- 800.0 


201.5 ± 


3.0 


513.0 ± 19.9 


22.4 ± 20.0 


050730 


300.0 - 


- 600.0 


245.6 ± 


2.4 


434.4 ± 16.3 


105.1 ± 16.5 


050730 


600.0 - 


- 800.0 


165.4 ± 


1.7 


672.2 ± 22.7 


54.1 ± 22.8 


051117A 


800.0 - 


1250.0 


365.0 ± 


1.7 


997.1 ± 50.0 


158.9 ± 50.1 


051117A 


1250.0 - 


1725.0 


389.3 ± 


1.4 


1328.3 ± 33.8 


41.5 ± 33.8 


060124 


300.0 - 


- 650.0 


221.7 ± 


1.6 


563 8 ± 7 9 


165.7 ± 7.7 


060124 


650.0 - 


- 900.0 


158.4 ± 


1.4 


694.9 ± 6.0 


28.5 ± 6.0 


060204B 


100.0 - 


- 270.0 


103.3 ± 


6.4 


118.6 ± 3.5 


13.2 ± 3.6 


060204B 


270.0 - 


- 450.0 


91.9 ± 


3.9 


332.3 ± 9.8 


33.4 ± 10.1 


060210 


165.0 - 


- 300.0 


98.5 ± 


1.1 


207.6 ± 6.5 


35.4 ± 6.4 


060210 


350.0 - 


- 450.0 


79.2 ± 


0.9 


369.9 ± 4.4 


10.4 ± 4.4 


060418 


83.0 - 


- 110.0 


22.7 ± 


0.4 


87.8 ± 1.2 


3.4 ± 1.2 


060418 


122.0 - 


- 200.0 


58.0 ± 


0.6 


130.3 ± 0.8 


4.8 ± 0.8 


060607A 


93.0 - 


- 130.0 


31.1 ± 


0.3 


99.0 ± 3.2 


3.7 ± 3.2 


060607A 


220.0 - 


- 400.0 


138.8 ± 


1.9 


265.3 ± 11.6 


35.5 ± 11.7 


060714 


100.0 - 


- 125.0 


15.9 ± 


0.2 


114.1 ± 0.6 


5.8 ± 0.6 


060714 


125.0 - 


- 160.0 


27.9 ± 


0.3 


132.5 ± 3.5 


5.2 ± 3.5 


060714 


160.0 - 


- 230.0 


48.3 ± 


1.2 


178.4 ± 1.7 


14.9 ± 1.6 


060904A 


250.0 - 


- 600.0 


219.6 ± 12.8 


288.9 ± 19.5 


21.0 ± 19.6 


060904A 


600.0 - 


1000.0 


314.2 ± 6.9 


678.5 ± 7.7 


33.3 ± 7.7 
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